Theory of Microwave- Assisted Supercurrent in Diffusive SNS Junctions 
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The observation of very large microwave-enhanced critical currents in superconductor-normal 
metal-superconductor (SNS) junctions at temperatures well below the critical temperature of the 
electrodes has remained without a satisfactory theoretical explanation for more than three decades. 
Here we present a theory of the supercurrent in diffusive SNS junctions under microwave irradiation 
based on the quasiclassical Green's function formalism. We show that the enhancement of the 
critical current is due to the energy redistribution of the quasiparticles in the normal wire induced 
by the electromagnetic field. The theory provides predictions across a wide range of temperatures, 
frequencies, and radiation powers, both for the critical current and the current-phase relationship. 



It was predicted by Eliashberg already in 1970 pQ that 
the condensation energy of a superconducting thin film 
can be increased by irradiating the film with microwaves. 
Within the framework of his theory, one can explain the 
microwave-induced increase of the critical current of su- 
perconducting bridges for temperatures very close to the 
critical temperature [2H1] > which is known as the Dayem- 
Wyatt effect. However, Eliashberg's mechanism fails to 
explain a related effect in diffusive SNS junctions. Several 
experiments have shown that upon irradiation the crit- 
ical current can be enhanced by up to several orders of 
magnitude, even at temperatures well below the critical 
temperature of the superconducting electrodes 0(5]. Ad- 
ditionally, these experiments show that the critical cur- 
rent is a nonlinear function of the radiation power, which 
existing linear response theories [8] cannot explain. 

There is now renewed interest in this problem, trig- 
gered by recent experiments. Fuechsle et al. [9j mea- 
sured the current-phase relationship under microwave ir- 
radiation, and reported that the current is progressively 
suppressed at phase differences close to 7r as the radi- 
ation amplitude increases. Moreover, Chiodi et al. [10] 
observed that critical current is enhanced when the mi- 
crowave frequency is larger than the inverse diffusion time 
in the normal metal. 

To understand the microwave-assisted supercurrent in 
diffusive SNS junctions, we develop a microscopic the- 
ory based on the quasiclassical Keldysh-Usadel approach, 
which takes into account the nonlinear effects of the mi- 
crowave irradiation. Our theory provides a quantitative 
description for a wide range of values of the temperature, 
microwave power, frequency, and the strength of inelastic 
scattering. In particular, we show that the large enhance- 
ment of the critical current originates from the presence 
of a minigap, E g , in the density of states of the normal 
wire. This minigap blocks some of the transitions caused 



by the microwave radiation, which results in a redistribu- 
tion of quasiparticles, enhancing the supercurrent when 
the temperature T is comparable or larger than Eg/hrs- 
We also show that the nonequilibrium distribution in 
the normal wire leads to a highly non-sinusoidal current- 
phase relationship, in a good agreement with Ref. [9]. 

We consider a diffusive normal metal (N) of length L 
connecting two bulk superconductors with energy gap A 
(see inset of Fig. [ljb)). In the absence of microwaves, 
superconducting pair correlations leak into the normal 
metal modifying its properties. For instance, the local 
density of states (DOS) is modulated [TT] and a super- 
current can flow through the normal metal [12]. The DOS 
exhibits a minigap E g (ip), see Fig. [lja), which depends 
on the superconducting phase difference ip [13 . For ideal 
interfaces, which we consider hereafter, E g (0) 3.12Et, 
where Et = hD/L 2 is the Thouless energy and D is the 
diffusion constant, whereas E g (ir) = 0. 

We model the microwave radiation by an oscillating 
electric field, E(i), described by a time-dependent vec- 
tor potential A(t) — Aq cos(woi), where Aq points along 
the axis of the junction. We neglect screening, and as- 
sume that the field is position independent []3] . We also 
neglect the effect of the radiation inside the superconduc- 
tors, which is justified for frequencies smaller than A/Ti 
or when the electrodes are thick compared to the size 
of the junction. To evaluate the physical observables, 
we use the quasiclassical theory of superconductivity for 
diffusive systems [T5J HE]- It is formulated in terms of 
momentum averaged Green functions G(R,t,t') which 
depend on position R and two time arguments. These 
propagators are 4x4 matrices in Keldysh/Nambu space: 

Here, G R ' A ' K are the retarded, advanced and Keldysh 
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FIG. 1: (Color online) (a) Local density of states (DOS) 
in the middle of the normal wire for A = IOOEt and dif- 
ferent values of the phase difference ifi/n, in the absence of 
microwaves, (b) Absorption rate 77+ for a high frequency 
Hujo/Et = 8 and (p — n/2, s — 0.125. Thin line shows the ap- 
proximation from Eq. Q. Inset: Schematic representation of 
the SNS junction, (c) Correction 6f = / — fo to the electron 
distribution function vs. energy at two different temperatures 
for tp = 7r/2, Tiuio/Et = 4, and s = 0.125. Solid lines cor- 
respond to the exact numerical results and the dashed lines 
to the approximation in Eq. p}. The thin black line shows 
the spectral supercurrent js(E) in the absence of microwaves, 
(d) The same as in (c) for ks/Er = 15 and <p = 0.87T. 

components, respectively. The Green functions acquire 
the BCS value Go{t — t',±ip/2) inside the superconduc- 
tors, and in the normal metal fulfill the Usadel equation 

HDVo (Go VoG) = [-ief 3 + ia,G} , (2) 

where o denotes the time convolution (X o Y)(t,t') — 
dtiX(t, tx)Y(ti,t'), V the gauge-invariant gradient 

V o X — VX — {[eAr^/h, X] which involves the vector 
potential A(t,t') = A(t)S(t - t'), and e(t, t') = ihd t 6(t - 
t'). The self-energy a describes inelastic interactions in 
the wire, and the Green functions are normalized as (Go 
G)(t,t') = S(t-t'). 

Because Andreev reflection blocks sub-gap heat trans- 
port out of the junction, inelastic interactions play an 
important role in balancing the effect of microwaves. We 
describe these interactions, for example due to phonons, 
within the relaxation time approximation. In this ap- 
proximation, the interaction strength is characterized by 
a constant energy (scattering rate) T [T7] . The microwave 
coupling, in turn, introduces the energy scale Ea = 
e 2 DAl/h. One can show that the ratio s 2 = Ea/Et 
determines the change in the spectral quantities due to 
the microwaves, while the ratio a — Ea/Y controls the 
corresponding change in the electron distribution. We 



note that s = eVo/hujQ, where Vq is the amplitude of the 
oscillating voltage across the junction. 

In order to solve Eq. ^ , we follow Ref. [18] and Fourier 
transform the Green functions to energy space. Due to 
the time dependence of the vector potential, the Usadel 
equation in energy space admits a solution of the type 
G(R, E, E') = J2 m Go,™(R, E)6(E — E' + mhuj ). With 
this Ansatz the Usadel equation becomes a set of cou- 
pled differential equations for the Fourier components 
G n>m (E) = G(E + nhujQ, E + mTkuo). For arbitrary radi- 
ation power we solve these equations numerically using a 
Jacobian-free Newton-Krylov method [T^l [2D] . From the 
solution of G, we can compute all physical observables. 

We now analyze the linear response regime (s 2 , a <C 1). 
In this limit, we can derive the kinetic equation for the 
time-average of the distribution function, /, by keeping 
terms up to second order in Aq in the Keldysh compo- 
nent of the Usadel equation ([2|. Because of Andreev 
reflection, when relaxation processes are slower than the 
diffusion inside the junction, / is in our gauge constant 
throughout the normal wire. Consequently, the kinetic 
equation reduces to an equality of the electron-phonon 
and microwave collision integrals, / -ph = Ij, averaged 
over the junction volume f2. The microwave collision in- 
tegral resembles Joule heating and it is proportional to 
a time-averaged product of electric field and current (at 

energy E), Z 7 = f^E(t) ■ Ttt 3 j k (E + hu /2, t) - (E ^ 

E - hw ), where j K = G R o V o G K + G K o V o G A . 
Using this result, the kinetic equation for the correction 
8f = / — fo, where fo is the Fermi function, becomes 

T(p)Sf = V-(E + ftwo)/+(l - fo) - V+(E)fo(l - /+) 

+ r,+ (E - &jo)/_(1 - fo) - V-(E)f (l - /_) . (3) 

Here, (p) is the spatially averaged density of states inside 
the junction and f± = fo(E ± hu)o). The emission 
and absorption (rj + ) rates are defined as T]+(E) = rj-(E+ 
hcu ) = - e -^lm Tr f 3 (E))/(f (E + %u )- fo (E)). 

For frequencies hu>o < 2E g (ip), one can neglect the AC 
components of the retarded/advanced functions, so that 

V+ « ^{p oP+ + Re{(/c? + ft)(f? + /+*)}) • (4) 

This reduces to the original linear response result by 
Eliashberg in the case of a bulk superconducting film 
PP. One can now see that the minigap in p and f R 
blocks some of the radiation-induced transitions (see 
Fig. [TKb)) . Thus, if the temperature is sufficiently high 
{k^T> E g ((p)), an excess of quasiparticles accumulates 
below the minigap, and their number is depleted above it. 
This cooling effect is illustrated in Fig.Qc), where the re- 
sult of Q is compared to the exact numerical result. As 
one can see, Eq. Q reproduces the main features of the 
exact result well especially at ksT ^> E g (cf>). Note that 
despite the cooling at some energies, the Joule power ab- 
sorbed in the junction, P = v F f dE EQ{I^) = I(t)V(t) 
where vp is the normal state DOS, is positive. 
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FIG. 2: Correction to the critical current, normalized by the 
critical current in the absence of the field, as a function of 
temperature and inelastic rate for A = 100-Et and several 
frequencies. The field strength is s = 0.125 in all cases. The 
lines separate the region of parameters for which the critical 
current is enhanced from that in which it is reduced. 

However, Eq. Q does not describe correctly the be- 
havior of the distribution function when Tiloq > 2E g (p) , 
as shown in Fig. Wd). This means that Eq. Q always 
fails to describe the behavior close to ip = ir. In this 
limit, the radiation induces changes in the AC compo- 
nents of the retarded/advanced quantities that couple to 
the time-averaged distribution function, especially at en- 
ergies close to E = ±hui /2. Since the behavior of these 
components is determined by a complicated balance be- 
tween diffusion and AC excitation, an accurate descrip- 
tion of rj± in general requires a numerical calculation. 

In the limit Y -C Ex, the correction to the supercur- 
rent comes mainly from the change in the distribution 
function, and it can be written as 

f° dEj s (E)6f(E), (5) 

where S is the cross-section, er/v = ^vpD the normal- 
state conductivity of the wire and js (E) the equilibrium 
spectral supercurrent [21], which is plotted in Figs. [TJc) 
and [ljd) together with Sf. Based on this, the cooling 
effect described by Eq. Q is expected to manifest as 
an enhancement of the critical current for fee? 1 ^ E g (0). 
This is confirmed by the exact numerical calculations ob- 
tained in the low-amplitude regime, see Fig. [2] The effect 
increases up to frequency Huiq « 2E g (0), and at larger fre- 
quencies becomes more varying, due to the complicated 
energy dependence of js and /. On the other hand, as Y 
increases, the magnitude of Sf decreases, which together 
with the suppression of minigap and js reduces the cur- 
rent. The above is in qualitative agreement with existing 
experiments which concentrated on huj /E T < 10. 

For high power, the magnitude of the critical current 
eventually decreases as can be seen in Fig. [3] This oc- 
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FIG. 3: (Color online) Critical current (normalized by the 
current without AC field) versus radiation amplitude s for a 
wire with A/Et = 100. (a) For different temperatures at 
TujJq/Et = 4 and Y/Et = 0.05. (b) For different frequencies 
and k B T/E T = 15. 

curs as large-amplitude oscillations of phase average the 
density of states, which results in a suppression of the 
coherence and in the subsequent closing of the minigap 
(see Fig. ga)). As a consequence, the cooling effect is 
suppressed, and microwaves mainly heat the electrons 
in the same way as in the normal state, which reduces 
the current. In the relaxation time approximation the 
temperature is for high field strength given by T* 
[P/(2iy F nk 2 B Y Q )} 1 ^ (provided T < T* < A/k B and as- 
suming Y(T) = YY T 3 [22 ), where P = a N nAlw%/2 is 
the average Joule power dissipated in the junction. 

The critical current also exhibits oscillations when ra- 
diation amplitude increases, see Fig.ga), similar to those 
already seen in the early experiments [6]. For short 
junctions (A < Et', not plotted), we find that these os- 
cillations match reasonably well with the usual Bessel os- 
cillations in Josephson tunnel junctions, i.e., / cx Jo (2s), 
but in the long-junction limit the similarity is only quali- 
tative. Locations of the dips in the I c (s) relation are not 
strongly dependent on the temperature, but depend on 
the radiation frequency, as shown in Fig. gb) . 

The microwave irradiation alters the current-phase re- 
lationship, enhancing the current at ip < | and suppress- 
ing it or even changing its sign at ip > | (see Fig. gb)). 
The behavior near tp = tt comes from two sources: the 
cooling disappears as the minigap closes, and the features 
peculiar to the dissipative AC response of SNS junctions 
not contained in Eq. Q become increasingly important. 
One can for example see in Fig. [ljd) that at ip/n = 0.8 
the peaks at E = ±hujo/2 give a significant negative con- 
tribution to the current. 
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current depends on the field strength non-monotonically, 
and predicted the modification of the current-phase re- 
lation. Moreover, our results pave the way for filling 
some remaining gaps in the understanding of SNS junc- 
tion physics such as the finite-voltage Shapiro steps [H] 
or the role of phase fluctuations providing the "intrinsic 
shunting" [25] . 
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FIG. 4: (a) Distribution function (solid black) and density 
of states (green) for large amplitude s = 2, and Hljo/Et = 4, 
A/E T = 100, k B T/E T = 10, T/E T = 0.05, and ip = tt/2. 
Microwaves cause heating from ksT / Et = 10 (dotted) to 
fc s T « hujo^ Ea/4F (blue), (b) Current-phase relation nor- 
malized to equilibrium critical current at ksT/Er = 15, 1 
(top to bottom) and s = 0.125 (solid) and 0.25 (dashed), 
for hcoo/Er = 4, A/E T = 100, F/E T = 0.05. (c) Current- 
phase relation for different amplitudes s = 0, 0.2, 0.3, 0.5, 0.75 
(solid, top to bottom) at ksT '/ ' Et = 10 and Hljo/Et = 1.2. 
Relaxation rate is chosen as T/Et = 0.2(fesT/10i5T) 3 , and 
A/Et = 58. Experimental data from Ref. [91 is shown as 
dots, (d) As in (c), for s = 0.3, Hujo/Et ~ 2, and tempera- 
tures k B T/E T = 8, 9.5, 10, 11, 12.5 (top to bottom). 

To compare with the results of Ref. [5] , we compute the 
current-phase relationship using the parameters of the 
experiment (T, Et,ojq). We have two free parameters: 
T/Et, which we assume large enough to suppress the en- 
hancement of the critical current, and the amplitude of 
the AC bias, which we fix by assuming that s — 0.5 cor- 
responds to the externally applied power level 28 dBm 
at Hujo/Et = 1.2. The power dependence, see Fig. |4|c), 
reproduces the main experimental features: (i) with in- 
creasing power the maximum supercurrent is reached at 
'/'max < tt/2, (ii) the supercurrent is strongly suppressed 
for phases close to n, and (iii) for ip < (p max , the super- 
current is slightly enhanced compared to s = 0. On the 
other hand, as shown in Fig.[4]jd), the deviation from the 
sinusoidal form becomes slightly more pronounced as T 
increases, in a qualitative agreement with experiments. 
The difference to the experiment at high power or low 
temperatures may be due to nonlinear radiation coupling 
and the relaxation time approximation, respectively. 

In summary, we have presented a general theory for de- 
scribing the effects of radiation on the properties of diffu- 
sive SNS junctions, which explains a wide range of exper- 
imental observations. We have clarified the mechanism 
of stimulated superconductivity, shown how the super- 
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